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Abstract 

A model calculation for pion production in nucleon-nucleon collisions is presented. 
Direct production, pion rescattering and contributions from pair diagrams are taken 
into account. The amplitudes for the elementary processes are based on well es¬ 
tablished microscopic models of nucleon-nucleon and pion-nucleon scattering. The 
A(1232) is included explicitly and is found to play an important role even at en¬ 
ergies close to production threshold. A good overall agreement with existing data 
from the pion production threshold up to the Delta resonance region is achieved. 


Over the last decade or so a wealth of rather accurate data on pion production 
in nucleon-nucleon ( NN ) collisions near threshold has become available. This 
concerns reaction channels with a two-body final state such as pp —> dn + [1- 
3], as well as channels with three-body final states such as pp —* ppn° [4-6] or 
pp —> pnir + [8,9]. The same period of time also witnessed lively activitiy on the 
theoretical side. In this case, however, most of the efforts concentrated on the 
study of one particular reaction, namely the process pp —> pp7r°. Clearly this 
development was initiated by the unexpected failure of the first model studies 
[10,11] to reproduce the corresponding data [4], The underprediction of the 
pp —> ppir 0 total cross section by a factor of 5 turned out to be a challenging 
problem for theorists. 

These first investigations followed very closely the approach of Koltun and 
Reitan [12]. In this model it is assumed that the pions are produced either 
directly from the nucleon (Fig. la) or via pion rescattering (Fig. lb), where in 
the latter the ttN amplitude is approximated by the corresponding tiN s-wave 
scattering lengths. Subsequently two ’’new” mechanisms have been suggested 
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that, in principle, allow one to obtain quantitative agreement with the data 
on 7 r° production close to threshold. In one of them the pion is produced 
via an intermediate virtual antinucleon-nucleon state in conjunction with the 
exchange of a heavy meson (p, lu, a, ...)(HME), as depicted in Fig. lc. It 
was found in Refs. [13,14] that the contributions from a and u exchange can 
enhance the production cross section significantly enough to reproduce the 
data. In the other mechanism the static nN interaction (i.e., the on-shell nN 
amplitude at threshold) used in the Koltun-Reitan model for the rescatter- 
ing process is replaced by an off-shell nN amplitude. Physically, the latter 
is required anyway because the exchanged pion in the rescattering process is 
off-mass-shcll. Since the isoscalar part of the nN amplitude, which is respon¬ 
sible for rescattering in the 7r° production, is practically zero on-shell and at 
threshold clue to constraints of chiral symmetry, but can be fairly large once 
one goes off-shell, the use of an off-shell amplitude, given by a certain model, 
leads likewise to an appreciable enhancement of the production cross section 
[15-17]. 

Despite these apparent successes there is still a controversy over which of the 
two processes if any, is the physically correct one. The HME mechanism has 
been put into question in a recent paper by J. Adam et al. [18]. These authors 
work in the context of the relativistic Gross equation where contributions 
from intermediate antinucleon-nucleon states are generated automatically by 
the scattering equation and summed up to all orders. They see only moderate 
effects from these virtual intermediate states, resulting in contributions that 
are considerably smaller than the ones obtained in the perturbative treatment 
of Refs. [13,14], Likewise, results based on the off-shell nN amplitude have 
come into dispute. Calculations done in lowest order chiral perturbation theory 
[19-21] give rise to a nN amplitude that differs in sign from the one employed 
in Refs. [15-17] when extrapolated off-shell. Accordingly, the amplitude from 
the rescattering diagram also changes sign and then interferes destructively 
with the contributions from the other production mechanisms. Consequently, 
agreement with the data in this case can be no longer achieved. For a thorough 
discussion of this topic see also Ref. [22], Furthermore we want to call attention 
to a recent paper by Bernard et al. [23], where it is argued that heavy baryon 
chiral perturbation theory (on which the works [19-21] are based) can not be 
used in reactions with as large momentum transfers as they are typical for 
meson production in nucleon-nucleon collisions. 

One possible way to learn more about the role of the individual production 
mechanisms has been suggested in Ref. [17,22], It consists of a comparative 
study of all relevant pion production channels (pp —> ppn°, pp —> pnn + , 
pp —> dn + and pn —> ppn~) in a consistent approach (i.e., with the same 
model and same parameters). Moreover, higher partial waves should be taken 
into account. Most of the aforementioned investigations consider only the low¬ 
est partial waves in the outgoing channel, which means that the NN system 
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is in an S-wave state (or in the deuteron bound state) and the pion is like¬ 
wise in an s-wave relative to the nucleon pair. Such calculations permit only 
conclusions on the absolute magnitude of the production cross section near 
threshold. The inclusion of higher partial waves, in the NN as well as the nN 
sector, allows one to calculate predictions for differential cross sections and, 
in particular, spin-dependent observables. Therefore it is possible to exam¬ 
ine whether the considered production mechanisms lead to the proper onset 
of higher partial waves, as suggested by the data. Results of the latter kind 
are particularly interesting because they reflect the spin-dependence of the 
production processes and therefore should be very useful in discriminating 
between different mechanisms. 

In this letter we present a model calculation for the reactions pp —► ppn°, 
pp —> pnn + 1 pp —> dn + , and pn —* ppn~. It is an extension of our earlier study 
for s-wave pion production [16,17] to higher partial waves. Now all NN- and 
7 tN partial waves up to orbital angular momenta L = 2, and all states with 
relative orbital angular momentum / < 2 between the NN system and the pion 
are considered in the final state. Furthermore, the excitation of the A(1232) 
resonance is taken into account explicitly. Therefore we are able to calculate 
meaningful predictions for the different reaction channels of NN —» NNn 
from their threshold up to the A resonance region. 

The reaction NN —> NNtc is treated in a distorted wave born approximation, 
in the standard fashion. The actual calculations are carried out in momentum 
space. For the distortions in the initial and final NN states we employ the 
model CCF of Ref. [24], This potential has been derived from the full Bonn 
model [25] by means of the folded-diagram expansion. It is a coupled channel 
(NN, N A, AA) model that treats the nucleon and the A degrees of freedom on 
equal footing. Thus, the NN ^ N A T-matrices that enter in the evaluation of 
the pion production diagrams involving the A isobar (cf. Fig. 2) and the NN 
T-matrices that are used for the diagrams in Fig. 1 are consistent solutions 
of the same (coupled-channel) Lippmann-Schwinger-like equation. 

The nN —> nN T-matrix needed for the rescattering process is taken from a 
microscopic meson-exchange model developed by the Jiilich group [26]. This 
interaction model is based on the conventional (direct and crossed) pole di¬ 
agrams involving the nucleon and A isobar as well as t-channel meson ex¬ 
changes in the scalar (a) and vector (p) channel derived from correlated 2n- 
exchange. Note that in our model of the reaction NN —> NNn contributions 
where the pions are produced directly from the nucleon or A (cf. Figs.la and 
2a-c) are taken into account explicitly. Therefore, the corresponding nucleon 
and A pole terms have to be taken out of the nN T-matrix in order to avoid 
double counting. 

For the nNN and nNA coupling constants at the pion production vertices 


3 



we take the values /^ 7V7r /47r = 0.0778 [24] and f^An/^ 7r = 0-26 [27]. The form 
factors at these vertices are chosen to be soft (We use a monopole form with 
a cutoff mass A*. = 900 MeV) in line with recent QCD lattice calculations 
[28] and other information [29,30]. The width of the A isobar is taken into ac¬ 
count by using a complex A energy in the propagator. Specifically, we employ 
a parameterization of the width given by Kloet and Tjon in Ref. [31] which 
is energy- as well as momentum-dependent. We wish to point out, however, 
that the width influences the observables only for energies Tj \ a f, > 420 MeV. 
For lower energies the results obtained with and without A width are prac¬ 
tically the same. Note that for simplicity we have suppressed the three-body 
singularity that appears in the pion rescattering diagram for energies above 
threshold by fixing the corresponding three-body propagator to its threshold 
value. Thus, three particle unitarity is not fulfilled in our calculation. Ear¬ 
lier studies [32] have shown, however, that a large part of the imaginary part 
produced by the three particle cut of the pion exchange is cancelled by the 
imaginary part arising from the nucleon self energy contribution above pion 
threshold. Therefore we expect the effect of our approximation to be small, 
at least for energies close to threshold. Furthermore, a non-relativistic boost 
is applied for the 7rA-T-matrix. We expect this to be appropriate for the 
s-waves. The effect of this treatment for the higher partial waves in the 7rN- 
systern needs further study. We come back to this point when we discuss our 
results. 

Results for total cross sections in all experimentally accessible channels are 
shown in Fig. 3 as a function of r), the maximum momentum of the produced 
pions in units of the pion mass. Evidently the predictions of the model are in 
good overall agreement with the data over a wide energy range. We emphasize 
that the results for the channels pp —> pnn + and pp —> dn + do not involve any 
adjustable parameters and are, therefore, genuine predictions of our model. In 
the case of pp —> ppre 0 , however, the basic model (including direct production 
plus rescattering) yields only about 60% of the measured cross section. (Cor¬ 
responding results are indicated by the dash-dotted curve in Fig. 3.) Here we 
have added contributions from the HME mechanism, Fig. lc, and fixed their 
’’strength” so that we can reproduce the data in the near-threshold region 
(cf. Ref. [17]). Specifically, we have included contributions due to u exchange 
using the vertex parameters NN /^n = 10 and A uNN = 1-5 GeV (monopole 
form factor). We would like to emphasize, however, that we do not view our 
HME contribution as being due to a genuine process but rather as an effective 
parametrisation of short range mechanisms [33] not considered explicitly. 

Note that, compared to our earlier work [16,17], now both time orderings of 
the rescattering diagram (Fig. lb as well as 2d and e) are properly included 
(the one where the pion is emitted off one nucleon and interacts with the other 
before emission as well as the one where one nucleon emits two pions, one of 
which is absorbed on the other nucleon). This leads to a considerable enhance- 
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ment of the rescattering contribution so that now no additional contribution 
from the HME mechanism due to the o meson is needed. The effect of the 
HME contributions on the reaction channels pp —► dn + and pp pnn + is neg¬ 
ligible [17,34], so that the corresponding results remain practically unchanged. 
Therefore we do not show them separately in Fig. 3. 

Let us discuss the influence of the A resonance on the production cross sections 
close to threshold. Our model includes pion production from the A directly 
(Fig. 2 a-c) or via nN rescattering (Fig. 2 d,e). The latter clearly gives con¬ 
tributions to s-wave pion production. In this context we wish to point out 
that now - unlike the purely nucleonic case - charge-exchange rescattering 
is possible even in the reaction pp —► ppn 0 via a A ++ n intermediate state. 
However, it is less known that also direct pion production from the A gives a 
non-zero contribution at threshold. This can be easily seen from the standard 
reduction of the nN A vertex starting from the Lagrangian 

CnNA = —+ h.C. , ( 1 ) 

rru 


where ip, $ and ip^ denote the nucleon, n and A field operators, respectively, 
which leads to the following expression: 
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Here p (F A = \J + p 2 ) is the momentum (energy) of the incoming A, q 

(t J q = \jm\ + q 2 ) the momentum (energy) of the produced pion and S the 
spin transition operator. Evidently, even for vanishing pion momentum q, the 
term proportional to uj q /M A , which is the analog to the recoil term appearing 
in the iViV7r-vertex, survives. 

Results for the production cross sections without inclusion of the A isobar 
are indicated in Fig. 3 by the dashed lines. These curves are obtained by 
setting the tiNA coupling in the production operator to zero. Obviously, for 
the reaction pp —> ppn 0 the contributions involving the A lead to a decrease 
of the cross section in the near threshold region. This reduction (by about 
20 %) is entirely due to the direct production mechanism; the contribution 
from rescattering off the A is negligibly small. At first this is very surprising, 
especially because - as was mentioned before - the A ++ intermediate state 
allows also rescattering in the dominant nN isovector channel. However, a 
detailed inspection of our results reveals that the rescattering contribution is 
only small because of a strong cancellation between the diagrams Fig. 2d and 
e. Individually their magnitudes are quite significant. 
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Also for the reaction pp —> dir + we find that the direct production is the domi¬ 
nant contribution among the ones involving the A. It increases the production 
cross section in the threshold region by about 30 %. Thus, its contribution is 
partly responsible for the observed overestimation of the dir + cross section 
close to threshold. Still we should stress that the contributions from the A 
that we get in our model are moderate as compared to the ones reported by 
Niskanen [34], In his case the inclusion of the A leads to an increase of the dir + 
cross section close to threshold by almost a factor of 3. We believe that this 
difference is due a different treatment of the three particle propagator in the 
rescattering contribution involving the A-resonance. Niskanen fixes the energy 
of the exchanged pion by putting it on-shell [35]. This choice maximizes the 
contribution of the A [36]. 

Analyzing powers for the reactions pp —> d7T + , pp —> pm: + and pp —► ppn 0 
at some selected energies are shown in Fig. 4. Evidently the data for this 
polarization observable are nicely reproduced by our model. This means that 
the model predicts the correct onset of higher partial waves, especially of p- 
waves, and also the correct ratio of p-waves to s-waves. The dashed curve in 
Fig. 4 shows the results without contributions involving the A isobar. It is 
clear that the inclusion of the A isobar is essential for reproducing the data. 
It plays an important role even for energies very close to threshold (r/ < 0.25). 

Double polarization observables, namely the spin-dependent total cross sec¬ 
tion A o T lo to t = — (A xx + Ayy) and the spin correlation coefficient A xx — A yy , 
have recently been measured at IUCF for energies r/ > 0.56 [37]. In Fig. 5 
we present the predictions of our model for these observables. Obviously the 
description of these data is not that good. Especially in case of A xx — A yy 
our model overestimates the data by a factor of about two. Since the result 
for A xx — A yy without the A (dashed line) looks much better one could get 
the impression that our treatment of the A-resonance is incorrect. However, 
this observable is also very sensitive to the non-resonant part of nN p-wave 
rescattering. To illustrate this we show in Fig. 5 a calculation, where all con¬ 
tributions from nN p-wave rescattering are switched off (dotted line). Also 
in this case the description of the A xx — A yy data clearly improves whereas 
the results for A cr T /a tot and A y remain practically unchanged. Thus, it seems 
that A xx — A yy is a particularly interesting observable for learning more about 
the relation between the resonant and non-resonant contributions from nN 
p-wave rescattering. 

With regard to our model calculation we have already mentioned that we use 
an approximate description for the boost of the nN T-matrix. It is possible 
that this approximation leads to an overestimation of the contributions from 
TiN p-wave rescattering. Further studies in this direction are required. 

In summary, we have presented a model for pion production in nucleon- 
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nucleon collisions where the production operator is derived in a framework 
consistent with the interaction potentials that are used for generating the 
amplitudes in the elementary (NN and nN) processes. The NN interaction 
includes explicit coupling to the NA channel so that a consistent evaluation of 
pion production from NN and N A states can be done. The nN amplitude is 
taken from a meson-exchange model of nN scattering developed by the Jiilich 
group. A good overall description of cross section data for the reaction chan¬ 
nels pp — > ppn°, pp —> pnn + , pp —> dn + and pn —► ppn~ from the threshold up 
to the A resonance region is achieved. Quantitative agreement, not only with 
integrated cross sections but also with analyzing powers, is found over a wide 
energy range. Thereby the inclusion of nN rescattering as well as of the A 
degree of freedom plays an important role. We also demonstrated that polar¬ 
ization observables are a powerful tool to investigate details of the dynamics 
of the process. 
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Fig. 1. Pion production diagrams taken into account in our model - nucleonic con¬ 
tributions. 



d) e) 

Fig. 2. Pion production diagrams taken into account in our model - contributions 
involving the A- resonance. 
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*1 = q*(max/ m * 


Fig. 3. Total cross section for NN —> NNir in the different charge channels. The 
dash-dotted line shows the result of the coherent sum of the direct production and 
the rescattering. For the solid line heavy-meson-exchange contributions are added as 
described in the text. The dashed line shows the result without contributions involving 
the A isobar. 
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pp^>d;r + 



pp->pn7i + 




Fig. 4. The analyzing power A y . The solid line is the result of the full model, the 
dashed line is the result without the A-isobar. Upper panel: pp —> dir + at Tfab = 
290. 7, 330 and f25 MeV, respectively. Experimental data are from Refs. [38,39]. 
Middle panel: pp —» pmr + at Ti a b = 300, 320 and 330 MeV, respectively. Exper¬ 
imental data are from Refs. [9]. Lower panel: pp —> ppn 0 at TL a b = 310, f80 and 
530 Me V, respectively. Experimental data are from Ref. [7]. 
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Fig. 5. Predictions of our model for the spin-correlation functions of the reaction 
pp —> ppir °. The solid line is the result of the full model. The dashed line is the 
result without the A -isobar, whereas the dotted curve shows the results, when the 
contributions from non-resonant irN p-wave rescattering are switched off. The data 
are from Ref. [37]. 
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